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Dissociating  Gases  as  Natural  Convection  Coolants1 
DANIEL  E.  ROSNER2 


Introduction 


The  use  of  liquids  near  their  critical  point  has  been  investigated  by  F.  Schmidt 
[l]  for  cooling  turbine  blades  by  natural  convection.  This  method  takes  advantage 
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of  (i)  buoyancy  induced- convection  currents  due  to  the  very  large  accelerations 
associated  with  turbine  blade  rotation,  and  (ii)  the  augmented  thermal  expansion 
coefficient  and  thermal  conductivity  of  liquids  near  their  critical  point.  Now,  it  is 
also  known  that  the  thermal  dissociation  of  polyatomic  gases  greatly  increases 
their  expansion  coefficient  [Z]  and  thermal  conductivity  [3-6],  It  therefore  seems 
relevant  to  inquire  as  to  whether  equilibrium  dissociating  gases  might  not  offer 
some  advantages  in  natural  convection  cooling  applications.  3  For  example,  in  the 
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case  of  a  turbine  one  might  ask  whether  comparable  cooling  effects  can  be  achieved 
with  less  total  mass  of  coolant  in  the  blade,  thereby  reducing  centrifugal  stress 
levels.  Our  purpose  here  then  is  to  present  and  illustrate  a  method  to  examine  the 
possible  merits  of  equilibrium  dissociating  gases  for  this  general  class  of  natural 
convection  cooling  applications. 

Conventional  Natural  Convection  Cooling  Systems 

A  typical  physical  configuration  envisioned  here  is  shown  in  Fig.  1.  The  coolant 
is  confined  in  one  or  more  holes  drilled  (in  the  spanwise  direction)  in  each  turbine 
blade  and  closed  at  the  outer  end.  The  centrifugal  acceleration,  a,  for  such  blading 
is  of  the  order  of  107  to  108  cm/sec2  and  is  thus  some  104  to  10s  times  greater  than 
the  gravitational  acceleration,  g. 

In  the  presence  of  temperature  gradients  produced  by  contact  with  the  convectively 
heated  blade  material,  the  coolant  expands  and  experiences  a  large  body  force  directed 
toward  the  blade  root.  Heat  is  removed  at  this  location  by  a  separate  cooling  system 
and  the  heavier  (coolant)  medium  now  experiences  a  body  force  driving  it  centrally 
toward  the  blade  tip.  Being  closed  at  this  end,  the  fluid  is  then  heated  once  again, 
turns  the  corner  and  repeats  the  cycle.  In  the  thin  thermal  boundary  layer  limit4 
attainable  heat  transfer  coefficients  have  been  estimated  by  Schmidt  [l]  for  the  case  of 
water  as  the  coolant  by  invoking  a  correlation  formula  of  the  Nusselt  type  Nu  = 

Nu(Gr,  Pr),  developed  from  observations  of  heated  vertical  plates  in  air.  Here  Gr 
is  the  familiar  Grashof  number,  defined  by 

Gr  E  a  p  •  AT  •  x3/v2  (l) 

4i.  e. ,  non-interacting  boundary  layers,  small  curvature  effects. 
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where  a  is  the  effective  acceleration,  (3  is  the  thermal  expansion  coefficient  of  the 

fluid,  AT  is  the  imposed  temperature  difference,  x  is  the  distance  from  the  origin 

of  the  thermal  boundary  layer,  and  v  is  the  kinematic  viscosity.  As  pointed  out  by 

Ellerbrock  [8]  Grashof  numbers  in  the  holes  of  free  convection  turbines  based  on 

the  length  of  the  holes  (x  =  L)  are  in  the  approximate  range  1010  to  10ls.  Since 

transition  to  turbulence  is  ordinarily  associated  with  values  of  Gr  •  Pr  of  the  order 

of  109  (the  so-called  critical  Rayleigh  number,  Ra  .  )  this  indicates  that  in  these 

cnt 

applications  the  boundary  layers  must  be  turbulent.  Hence,  the  correlation  formula 
suggested  by  Eckert  and  Jackson  [9]  would  be  appropriate,  5  viz. 


+  0.  494  (Pr) 


Vs 


(2) 


Equation  (2)  is  expected  to  be  most  accurate  for  Prandtl  numbers  of  order  unity,  if 
the  property  variations  across  the  boundary  layer  are  small.  In  a  variable  property 
situation  reasonable  predictions  are  possible  if  one  introduces  all  property  values 
evaluated  at  a  reference  temperature  about  midway6  between  the  highest  and  lowest 
temperatures  appearing  in  the  system. 

If  in  this  class  of  applications  one  is  interested  in  ranking  fluids  as  to  their 
ability  to  transfer  heat  at  high  rates  in  a  system  of  prescribed  a,  L,  AT,  then  Eq. 

(2)  indicates  that  the  relevant  property  group  to  consider  is 

1  +  0.  494  ^  j  (3) 

5The  exponent  of  the  Prandtl  number  is  erroneously  given  as  ?/ 5  in  the  second  edition 
of  Ref.  10. 

^Actually,  the  work  of  Sparrow  quoted  in  Ref.  10  indicates  that  T+  should  be  chosen 
somewhat  closer  to  the  hot-side  temperature. 
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which  may  be  regarded  as  a  reference  heat  flux,  i.  e.  ,  the  heat  flux  pertaining  to  an 
acceleration  of  one  g,  unit  length  and  unit  AT,  in  any  consistent  set  of  units.7  It  is 
therefore  of  interest  to  compute  values  of  this  group  over  the  temperature  and  pres¬ 
sure  range  of  interest,  for  both  liquids  and  gases,  with  and  without  chemical  reac¬ 
tion.  Some  initial  results  of  this  type  will  discussed  below. 

Effect  of  Dissociation  on  the  Thermal  Conductivity  and  Expansion  Coefficient  of  Gases 

7, 

According  to  the  simple  kinetic  theory  of  ideal,  inert,  gases:  \«T',pccT  , 

3/  7 

v  x  t  '2p~1  and  a  =  T  '  2p"1,  so  that  for  such  gases  one  expects  the  approximate  behavior 

-  4/  n/ 

^ref  T~  /l0  (4) 

Therefore  q"  will  ordinarily  be  largest  at  high  pressures8  and  low  temperatures, 
ref 

being  almost  directly  proportional  to  the  density.  However,  for  gases  which  disso- 
ciate  (e.g.  ,  H2  —  2H,  N204  -  ZN02,  Cl2  -  2C1,  2HI  -  H2  +  I2,  (HF)6-  6HF) 

large  departures  from  this  simple  behavior  are  known  to  occur.  This  is  particularly 
true  in  the  case  of  the  heat  conductivity  \  and  expansion  coefficient  p.  By  compar¬ 
ison,  dissociation  introduces  relatively  small  changes  in  the  kinematic  viscosity  v, 
thermal  diffusivity,  a,  and,  hence  the  Prandtl  number  Pr  =  v/a.  A  rough  estimate  of 
the  thermochemical  augmentation  in  is  then  provided  by  the  simple  approximate 

Alternatively,  q"  may  be  regarded  as  the  average  heat  transfer  coefficient  for  a 
configuration  with  unit  length  and  unit  temperature  difference. 

8Stress  considerations  will,  of  course,  set  an  upper  limit  on  the  allowable  internal 
pressure.  Calculations  below  have  not  been  made  for  pressures  above  10z  atm. 


4 


TP-83 

result 


q  " 

ref,  dissoc 

q» 

ref,  inert 


(5) 


For  example,  consider  Figs.  2  and  3,  which  were  prepared  for  the  case  of  equilibrium 
dissociating  hydrogen  [5,  6].  At  p  =  0.  1  atm,  X  /  X  ~  9.  5  and  (3  /  (3  88  3.  1  so  that 

6(^  X  0^  X 

o 11  „  /a"  15  at  a  temperature  level  of  about  3250°K.  At  p  =  1  atm,  X  /  X  ~ 

1ref,  eq7  ^ref,  f  ^  r  eq f 

7.  5,  j3  /p  =2.8  so  that  q"  ,  /q1'  *  11  at  ?.  temperature  level  of  3750°K.  Thus, 

eq  f  ref,  eq  ref,  f 

dissociation  can  increase  the  reference  heat  flux  by  more  than  one  order  of  magnitude, 
with  the  greatest  augmentation  being  realized  at  low  pressures.  Of  course  the  temp¬ 
erature  level  at  which  these  large  effects  occur  (for  any  specified  pressure)  will  de¬ 
pend  on  the  strength  of  the  chemical  bond  being  broken.  Indeed  it  can  be  shown  that 
for  any  dissociating  diatomic  gas  the  maximum  enhancement  will  occur  at  the  temp¬ 
erature  at  which  the  degree  of  dissociation,  a  ,  is  near  0.  50.  In  terms  of  the 

eq 

equilibrium  constant  K^(T),  this  will  occur  at  the  "optimum"  temperature  T^  ^ 

defined  implicitly  by  K  (T  )  ~  (4/3)  p.  Conversely,  if  one  desires  operation  at 

P  opt 

some  prescribed  temperature  level  T  then  the  maximum  advantage  of  dissociation 

will  be  realized  at  the  pressure  level  p  ~  (3/4)  K  (T). 

opt  p 

In  the  light  of  this  background  information  and  the  detailed  data  given  in  Refs.  3, 

4,  11,  and  12,  one  can  now  compute  the  "figure  of  merit"  q"  for  various  coolants, 

ref 

over  a  range  of  pressures  and  temperatures. 

Relative  Effectiveness  of  Dissociating  Gases  as  Heat  Transfer  Media  in  the 
Turbulent  Natural  Convection  Regime 

Values  of  defined  by  Eq.  (3)  are  shown  plotted  in  Fig.  4  over  the  temperature 

o  o 

range  300  K  <  T  <  5000  K  for  the  particular  dissociating  gases  H2,  N204  and  air. 
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For  both  H2  and  N204  the  expansion  coefficient  was  computed  from  chemical  contri¬ 
bution,  C  ,  ,  to  the  molar  heat  capacity,  with  the  help  of  the  relation  [  5l 

p,  chem 


P 


T 


(Cp.  chen/R)  ] 

(aiT7rt)  J 


(6) 


Relevant  data  was  obtained  from  Refs.  11,  3,  4,  and  12,  respectively.  For  dis¬ 
sociating  air,  use  was  made  of  the  compressibility,  Z,  tabulated  by  Hansen  [12], 
together  with  the  relation 


P  = 


9  In  Z 
9  In  T 


const J 


(7) 


derived  from  the  definitions  of  (3  and  Z.  The  trends  anticipated  in  the  discussion  of 
the  previous  section  are  borne  out  (cf.  Fig.  4).  In  particular,  it  is  interesting  to 
note  that  at  p  =  1  atm  and  T  =  300°K,  dissociating  nitrogen  tetroxide  outperforms 
either  air  or  hydrogen  by  more  than  one  order  of  magnitude  on  a  heat  flux  basis 
Also,  as  expected,  both  hydrogen  and  air  exhibit  their  maximum  thermochemical 
contribution  at  much  higher  temperature  levels,  reflecting  the  higher  strengths  of 
the  bonds  being  broken  [H2,  (103  Kcal/mole),  02  (118  Kcal/mole)  and  Nz  (225  Kcal/ 
mole)  as  compared  with  N204  (  ~  14  Kcal/mole)  ].  While  data  for  N204  at  higher 
pressures  was  not  readily  available,  it  appears  that  N204  could  be  used  to  obtain 
values  of  q"  ,  larger  than  10”3  cal(cm)'  2(sec)~5  at  pressures  less  than  102  atm. 
However,  it  is  also  clear  from  Fig.  4  that  liquids  outperform  pressurized  dis¬ 
sociating  gases  on  a  heat  flux  basis  (cf.  data  [10]  for  saturated  Hg,  NH3  and  HzO  ). 
However,  additional  data  is  required  at  higher  temperature  levels,  and  the  role 

of  stress  limitations  in  properly  assessing  the  relative  merits  of  liquid  coolants 
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in  a  particular  application  remain  to  be  investigated.  It  would  also  be  of  interest 

to  consider  saturated  vapor/liquid  mixtures  in  this  regard,  with  some  attention 

paid  to  cases  in  which  the  vapor  dissociates  [  1 3 ].  Finally,  additional  study  of 

phenomena  near  the  critical  point  is  necessary. 

Required  Values  of  the  Reference  Heat  Flux 

An  estimate  of  the  value  of  q11  required  in  any  particular  case  can  be  obtained 

on  the  basis  of  the  following  illustrative  heat  balance  considerations.  Suppose  it  is 

desired  to  have  a  turbine  blade  operate  at  a  temperature  AT  below  the  turbine 

inlet  gas  total  temperature.  Then,  if  h  is  the  mean  heat  transfer  coefficient  for 

g 

the  outer  area,  A  of  the  entire  blade,  the  total  forced  convective  heat  transfer  rate  will 
o 

be 
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For  example,  consider  a  case  for  which  L  =  5  cm,  h  =  10”1  cal(cm)"  2(sec)-1(°K)-1, 

8 

A  /A.  =  2,  and  a/g  =  105.  Then,  if  one  desires  a  cooling  effect  of  A  T  ~  300°K  and 

o  1  g 

one  can-maintain  the  blade  root  at  a  temperature  such  that  AT  ~  AT  one  finds  that  q"  * 

g  ref 

mustbeofthe  order  of  lO^ca^cm^sec)'1.  As  sugg.e  sted  by  Fig.  4,  values  of  this  order 
of  magnitude  can  be  achievedwith  dis  sociating  gases  at  internal  pressures  below  102atm. 
Conclusions 

On  the  basis  of  the  considerations  outlined  above  the  following  conclusions 
have  been  drawn: 

a.  Molecular  dissociation  increases  the  effectiveness  of  gases  as  natural  con¬ 
vection  heat  transfer  media  largely  by  virtue  of  the  attendent  increase  in  thermal 
expansion  coefficient  (3  and  thermal  conductivity  X.with  the  latter  effect  being  the 
more  important  (cf.  Eq.  5).  Consequently,  in  temperature  limited  systems  an  optimum 
pressure  level,  strongly  dependent  on  the  strength  of  the  bond  being  broken,  will 

exist  for  maximum  coolant  effectiveness. 

b.  Coolants  in  the  dissociation  region  can  be  ranked  in  accord  with  a  figure 
of  merit  derived  from  the  thermophysical  property  groupings  which  enter  the 
expression  for  the  relevent  convective  heat  transfer  coefficient  (cf.  Eq.  3). 

c.  Fully  enclosed,  pressurized  dissociating  gases  may  be  capable  of  trans¬ 
ferring  heat  away  from  turbomachinery  components  at  rates  sufficient  to  enable 
the  use  of-  higher  inlet  total  temperatures.  Even  at  1  atm,  lg  and  temperatures 
not  far  from  300°K,  easily  dissociated  gaseous  coolants  can  be  found  which  are 
more  than  ten  times  as  effective  as  air. 

d.  A  simple  ranking  technique  of  the  type  presented  can  be  used  to  show  when 

dissociating  gases  are  competitive  for  a  given  natural  convection  cooling  application 

and, hence,  when  more  accurate  boundary  layer  calculations  or  experiments  are 
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required.  As  already  indicated,  very  few  accurate  natural  convection  solutions 
for  dissociating  gases  have  been  reported  in  the  open  literature  to  date. 
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COMPARISON  OF  LIQUIDS  WITH  DISSOCIATING  GASES 


